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SUMMARY 


A joint industry-government flight test in 1996 at Crows Landing, CA demonstrated that 
high deceleration rates can significantly benefit noise abatement approach procedures for 
the Sikorsky S-76 [1]. Deceleration increases aircraft angle-of-attack, thereby decreasing 
blade vortex interaction (BVI) noise emissions for the S-76 at descent rates exceeding 
600 fpm. It is believed that deceleration effects are likely to play a beneficial role for 
other rotorcraft. However, testing with large microphone arrays such as at Crows Landing 
is both impractical and costly. Therefore the capability of predicting noise levels for 
decelerating approach procedures is needed to support future noise abatement flight 
procedure development. 

The current effort addresses development of a noise modeling capability for decelerating 
approaches. The resulting technique employs discretization of the descent trajectory as 
multiple steady state segments for input to CAMRAD.Modl to predict rotor states for 
acoustic analysis. Deceleration is included by modifying the CAMRAD.Modl free flight 
trim options to allow trim to the specified acceleration components. 

The Sikorsky effort included the following subtasks: 

1. Develop and document a CAMRAD.Modl input deck for the S-76B main rotor. 
Also, develop an aerodynamic coefficient look-up table input file for the S-76B 
fuselage including control surface contributions. 

2. Exercise this input deck for several sample cases that match fixed operating 
conditions flown during the 1996 Crows Landing flight test. Correlate the 
predicted rotor trim state for each sample case with the measured rotor data. 
Assess accuracy of input decks through a combination of wind tunnel mode (rotor 
only) and free flight (complete aircraft) trim cases. 

3. Develop a standalone kinematic trajectory analysis tool that will define the 
trajectory as a one-parameter function of time or distance. Constant glideslope 
descents or constant rate-of-descent flight conditions will be included in the 
analysis. Fixed flight parameters will be input while the free parameters will be 
calculated. To incorporate unsteady flight effects, the flight path trajectory will 
then be discretized using the distance or time parameter into a sequence of flight 
conditions to be run using CAMRAD.Modl. At each discrete point in the 
trajectory, the inputs for CAMRAD.Modl will be extracted and written to a script 
file in namelist form. These namelists correspond to the .scr file for 
CAMRAD.Modl. From the trajectory analysis will come namelist inputs such as 
altitude, flight-path angle, vehicle velocity vector, and vehicle acceleration vector. 
The outputs will be in a form such that CAMRAD.Modl will be run with the input 
variable NCASES equal to the number of discretized points needed to accurately 
define the trajectory. Therefore, only one CAMRAD.Modl job will be submitted 
for a complete trajectory analysis. This tool can become another component of the 
TRAC analysis allowing a standalone capability for flight path modeling. 
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4. Using the results from step 3 above, the CAMRAD.Modl code will be used to 
predict the vehicle trim state and rotor loads for each discretized point along the 
flight trajectory. The CAMRAD.Modl free flight trim options will be modified to 
allow trim to the specified acceleration components. These components are input 
from the trajectory analysis code. A ‘low frequency’ approximation will be used in 
which the pitch acceleration is assumed zero. 

5. The CAMRAD.Modl input deck and an informal task report will be submitted to 
NASA. 

The CAMRAD.Modl input deck was developed by modifying an input deck provided by 
NASA Langley. Owing to the fact that a leased S-76B aircraft was used for the Crows 
Landing flight test, the comprehensive aircraft state data typically available from an 
experimental aircraft was not available for this effort. In addition, the available aircraft 
state data evidenced significant fluctuations for the level flight cases initially used to 
exercise the CAMRAD.Modl input deck. Hence additional correlations against other 
Sikorsky S-76 flight test data and predictions were used to evaluate the CAMRAD.Modl 
input deck. The predicted levels were provided by GEN HEL, a Sikorsky handling 
qualities model [2], 

After showing good correlations between the CAMRAD.Modl results and the flight 
test/GEN HEL data, the standalone kinematic trajectory analysis tool developed in 
Subtask 3 was implemented with the CAMRAD.Modl model to evaluate vehicle trim 
state and rotor loads along the flight trajectory for decelerating approaches. Two test 
cases were used to assess the results for a decelerating approach at constant rate-of- 
descent (ROD) and a decelerating approach at constant glide slope. The resulting 
comparisons with test data showed good agreement between the predictions and the flight 
test data, in particular for angle-of-attack. 

Comparisons of trim conditions for rotor alone and pitch and sideslip angle comparisons 
for full aircraft cases are presented in this report. A description of the standalone 
kinematic trajectory analysis tool and interface to CAMRAD.Modl is included. The 
flowchart for approach trajectory acoustic analysis is included here to illustrate the logic 
flow of the total acoustic analysis. 

Problem Statement 


The objective of this task was to develop the analytical tools required to model noise 
abatement approach flight paths flown by Sikorsky’s S-76B helicopter in a joint industry- 
government flight test at Crow’s Landing, CA in 1996. The completed tools will be used 
by NASA researchers to assess current capability in predicting noise footprints on the 
ground. Rotor trim state predictions were compared to wind tunnel rotor alone test data. 
Vehicle trim state predictions were compared to Sikorsky flight test data and Crow's 
Landing data. Flight trajectory cases consisting of decelerating approaches at constant 
rate-of-descent and a constant glideslope were used to predict free flight vehicle trim and 
compared with flight test data from Crow’s Landing. For this purpose a standalone 
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kinematic trajectory analysis tool was developed and implemented to interface with 
CAMRAD.Modl. In addition, vehicle acceleration was included in the CAMRAD.Modl 
trim controller for each segment of the quasi-steady analysis. 

A look-up table of aerodynamic forces and moments for an S-76B fuselage as functions 
of angle-of-attack and sideslip was synthesized from one-fifth scale model wind tunnel 
test data [3]. 

The following sections describe efforts and results involved in each subtask of this task. 
Listings of required FORTRAN routines are included. All comparisons with test data are 
included in this report. 

Information on CAMRAD.Modl can be obtained from [4], 



CAMRAD.Modl input files for the S-76B were obtained from NASA Langley. These 
files were compared with existing S-76B analytic rotor models at Sikorsky and updated 
accordingly. Rotor trim states obtained from rotor alone cases were compared with wind 
tunnel test data [5]. 

Fuselage aerodynamic loads and moments for fuselage and horizontal and vertical tails 
were obtained from a one fifth model scale wind tunnel test carried out at the UTRC main 
tunnel [3], The data were formulated as a look up table of loads and moments as 
functions of fuselage angle of attack and sideslip angle. The switch WBTTAB in the 
namelist NLBODY was turned on in CAMRAD.Modl according to the description in 
Reference 2. This information is referenced to the location specified by the parameters 
FSWB, BLWB and WLWB in namelist NLBODY. Higher order effects of the rotor on 
the body angle-of-attack, which is a function of the forward flight speed, were ignored in 
determining these airloads. In addition, the effect of the rotor on the tail angle-of-attack 
as a function of forward speed was also ignored. The data specifying the variation of 
vertical tail sideflow angle with aircraft yaw angle appeared to be unusually high as 
obtained from the wind tunnel data [3]. Therefore the vertical tail angle-of-attack was 
made equivalent to the aircraft yaw angle. This assumption didn’t seem to adversely 
affect the predictions obtained from CAMRAD.Modl . The aerodynamic loads and 
moments are tabulated for the following variations of fuselage angle-of-attack and 
sideslip (in degrees): 

Fuselage pitch: -10.0, -7.8, -5.6, -3.3, -1.1, 1.1, 3.3, 5.6, 7.8, and 

10.0 

Fuselage yaw: -20.0, -15.0, -10.0, -7.8, -5.6, -3.3, -1.1, 1.1, 3.3, 

5.6,7.8,10.0,15.0,20.0 
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Current setup for nacelle angles, aileron and rudder deflections are still included and 
dummy variations are included for these. These make up 4480 data points in the input 
aerodynamic data file. 


RESULTS 


Subtask 2 Correlation of Rotor Trim State with Test Data 

Test cases with fixed operating conditions were compared for this subtask. Rotor alone 
case were compared with full scale isolated S-76 rotor 40 ft by 80 ft wind tunnel data 
from [5], Sample case for a collective sweep chosen with a forward flight speed of 100 
kts, 293 main rotor rpm, and a shaft angle of -4 degrees. This is equivalent to an advance 
ratio of 0.23. Figure 1 is a comparison of the collective obtained from wind tunnel data 
and CAMRAD.Modl predictions obtained as a function of CJa. Fairly good agreement 
was obtained. Lateral (0| C ) and longitudinal (0| S ) cyclic variations for zero flapping are 
shown in Figures 2 and 3, respectively. Comparison of the longitudinal cyclic is better 
than that of the lateral cyclic although both are within generally expected comparisons for 
cyclic pitch values. 

CAMRAD.Modl predictions were compared with Sikorsky flight test data for steady 
level flight conditions. Crow's Landing data for level flight were found to be highly 
oscillatory and not deemed satisfactory for comparison with CAMRAD.Modl for the 
level flight comparison only. The results from the comparison with the flight test and 
GEN HEL [2] data are given in Figures 4 to 8 for level flight speeds from 60 to 150 kts. 
Owing to the nature of the data, the values of the dependent variables have been removed 
in the figures. However, it is clear from these results that the relative comparisons 
between CAMRAD.Modl, flight test data and GEN HEL simulations are reasonably 
close. Figure 4 is the fuselage pitch angle and the CAMRAD.Modl values are close to 
the flight test data. Figure 5 shows the fuselage yaw values with forward flight speed. 
Flight test did not obtain a yaw angle. However the GEN HEL values seem to be close to 
the CAMRAD.Modl values for this case. Very good results were obtained for the main 
rotor collective as shown in Figure 6. Comparisons with the flight test data and GEN 
HEL are included in this figure. Lateral and longitudinal main rotor cyclic variation are 
shown in Figures 7 and 8 respectively. Longitudinal cyclic CAMRAD.Modl predictions 
compare well with flight test data as shown in Figure 8. Lateral cyclic CAMRAD.Modl 
predictions agree well with GEN HEL simulations. Perhaps this indicates that both 
numerical simulations have the same level of difficulty in predicting lateral cyclic. 

Subtask 3 Standalone Kinematic Trajectory Analysis Tool 

This subtask has been defined to develop a procedure for the numerical analysis of noise 
generated by a helicopter during an approach to landing. The unique description of an 
approach to landing requires a description of the trajectory taken by the helicopter from 
the start of the approach until landing. This description can come from several possible 
sources such as flight test or flight simulation software. In this task, the flight path 
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trajectories are classified by descriptors such as constant glideslope or constant rate-of- 
descent. Because these descriptors concern only the trajectory itself, and are not vehicle 
specific, it is possible to use a very simple kinematic analysis to describe the aircraft 
trajectory. If any descriptors had been aircraft specific, a full blown rotorcraft simulation 
like GEN HEL would have been necessary. 

The procedure will use the CAMRAD.Mod 1 code for the prediction of quasi-steady 
vehicle trim and rotor loads. Quasi-steady here means that CAMRAD.Mod 1 will be used 
to analyze a sequence of trajectory snapshots. Each snapshot can represents a constant 
time or constant downrange distance increment. At each snapshot, the instantaneous 
vehicle velocities, attitudes, and accelerations from the trajectory analysis will be 
converted to a CAMRAD.Mod 1 input case and written to a file that will become the 
CAMRAD.Mod 1 script file. The script file will contain the sequence of runs which will 
completely describe the flight path trajectory. Therefore, only one CAMRAD.Mod 1 job 
need be submitted for each trajectory analysis. The description of the kinematic trajectory 
analysis now follows. 

Analysis Description 

The analysis allows several options to describe the trajectory, the end point data, and the 
trajectory discretization method. These options are as follow: 

Trajectory Options: 

1 ) Constant glideslope descent 

2) Constant rate-of-descent 

End Point Specification Options: 

1) Specify initial point location 

2) Specify initial point acceleration 

Trajectory Discretization Options: 

1) Constant time increments 

2) Constant downrange distance increments 

One assumption used in the determination of the trajectories is that the acceleration is 
constant throughout the maneuver. Using this assumption and the following relations the 
trajectory can be uniquely defined. 

5(/) = JV(ryr 

V(/) = V : +at 

where S is the distance along the flight path from the initiation of the maneuver, V is the 
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velocity or the individual velocity components, a is the acceleration, and t is the time 
from start of the maneuver. The figure below shows a simple line drawing of the path for 
the constant glideslope trajectory. 
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There are several possible combinations of flight path trajectory description and initial 
condition specification. Appropriate combinations which are available in the current 
analysis include: 1 ) constant glide slope, initial point location, 2) constant glide slope, 
initial point acceleration, 3) constant rate-of-descent, initial point location, and 4) 
constant rate-of-descent, initial point acceleration. Descriptions of the calculations 
involved for each combination are now given. In each case, the solution involves 
identifying the time history of the aircraft position, velocity, and acceleration. 



In this case, the proper initial and end point conditions which are given include: 


£/» Zfi V/' 

Now, the total time for the descent maneuver, t„ can be solved for using the expression 
for the known total distance for the maneuver. This is written as 


/, t, 

s, = J vfr Wr = | (v, + a T )cl T = v,/, +T at; = v,7, +\{at, )/, 
0 0 


but s t is known from the following expression: 



7 



The acceleration may be solved for using the following equation for velocity: 


v, = v, +at, 

a = ^ 

t, 

at, = v, - v- 

Now, .v, may be equated to solve for the total time 
s , =y(v>, +V,)/, 


Using the above with the following relations: 


sin6 H , = — — — — , cos 0 pp =-— 
s ' s < 

v » = v ( . cos 0 pp , v.j = v, sin0 f/ ,,v (/ = v, cosO^v^ = v, sin0 F/ , 

the final expressions for the discretized trajectory may now be written: 

/ 

x (‘ ) = - v , + J ( v ’„ +a x T)=Xi+ vj + \ aj 2 

0 

z{t) = z i +v :i t+ia.r 

t 

v A t ) = v ii +\a x dT = v xl +a x t 

0 

v. (r) — v_- + o.t 
a K = acos0 Ir 
a_ — «sin0 f/ . 



A valid set of end point data is: 


■J’ V < i > V :i> V ,r V ;n a 


In this case, the initial velocity components are known and the initial starting point must 
be determined. This requires a simple manipulation of the terms above to determine the 
total trajectory length after which the start point coordinates may be computed. The 
expressions below are used in this computation: 
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V o 2 

'’w + v;. 

v j = V v * + v i 

a 

_ t i v , + v ) 

A ' 2 

Now the initial coordinates are computed geometrically: 

X, = S, COS0 f/> = -s, — 
v i 

z,=z, +s, sin 6 hP =z f -s, -2- 

The instantaneous values for the trajectory may now be evaluated using the same final 
expressions from Case 1 above. 

Case 3: Constant Rate-of-Descent, Initial Point Location 
A logical set of end point conditions for this case is: 

Xj * ’ <• j ’ C ’ 

In this situation, the total time for the maneuver can be directly calculated: 

V. 

>•„ = fFZi 

We need to compute the x acceleration to determine the final x velocity. 

x f =0 = x i +j(v xi +aj)dT 
0 

o 

t; 

x , = ~ V J, ~ — 

which gives for the x acceleration and final jc velocity. 
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~ 2 ( x , + V J , ) 

V xf = V .,i +«/» 

Now the time dependent trajectory information may be calculated as follows: 

-*(')= v, +vj + ±aj 2 
z(t) = z, + v : t 
v,(/) = v u +aj 
v. = const 

Case 4: Constant Rate-of-Descent, Initial Point Acceleration 
The inputs for this case are: 


The starting point coordinates need to be calculated here. Several terms can be directly 
determined: 


a x = a 

_ ( V v/- r ») 

Now the starting point coordinates are computed: 


All the terms required to compute the trajectory as in Case 3 are now available. 


Trajectory Discretization 


There are two options for discretizing the flight path trajectory as listed above. The 
constant time increments are easily calculated from: 


A t = 


(nt- 1) 
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where nt is the total number of discrete time steps. For the constant distance option, the 
discrete times at which snapshots of the trajectory are to be taken must be calculated. 
These times may be computed from the equation for trajectory distance as follows: 

t t 

s(t ) = | v{x)dr = | (v, + az)dr = Vjt +\cit 2 

0 0 

For constant distance, the individual distances are determined from: 


Av - 


s , 

(nt-l) 


Now, t may be calculated for a given distance as: 


\at~ + Vjt — .v(/ ) = 0 

t -v, ±V? + 2as{t) 

a 

In this case, to move forward in time, the positive root is chosen. 

Figure 9 is a flowchart showing the logic used to obtain the CAMRAD.Modl solutions 
for a specified flight trajectory. These results will be subsequently used in an acoustic 
analysis that was not the subject of this task effort. Trajectories are classified as constant 
rate-of-descent or constant glideslope. Constant rate-of-descent trajectories can be 
described with equal range increments or as equal time increments. Trajectories defined 
as constant glideslope are described with equal time or range increments. These 
specifications are illustrated in Figure 10. 

Sublask 4 Prediction of Vehicle Trim State and Rotor Loads along Trajectory 

In order to accurately predict the external noise for a helicopter during an approach to 
landing, the time history of the vehicle operating state is needed over the descent 
trajectory. As part of the current task, a separate trajectory analysis described previously 
is used to generate CAMRAD.Mod 1 input for discrete points along the trajectory. The 
vehicle operating state is then solved for at each point and the output becomes available 
for the external acoustic analysis. One important thing to note about the approach 
trajectories is the fact that the vehicle is undergoing a deceleration. Until now, 
CAMRAD.Modl had been configured to trim to zero net accelerations for the free flight 
trim options (OPTRIM=l-6). This was modified to allow trim to a specified acceleration 
in the vertical and/or horizontal directions. The modifications are meaningful only for the 
free flight trim mode. 
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The modifications were implemented, from a user standpoint, through the addition of two 
new inputs in the namelist NLTRIM. These are the accelerations in the inertial reference 
frame, AXTRIM for the trimmed acceleration in the x-direction and AZTRIM for the z- 
direction. They are activated by loading values for them into the input file. To return to 
the original form of CAMRAD.Modl, set them equal to zero. The units for the 

-) o 

accelerations are either ft/s' or m/s*, whichever is appropriate. 

The free flight trim equations in CAMRAD.Modl are written for a body-fixed frame of 
reference. The acceleration terms, which are more convenient to specify in the inertial 
frame from the trajectory analysis, must be transformed to the body reference frame. The 
transformation is described by the aircraft pitch and roll Euler angles. This 
transformation is already available in CAMRAD.Modl and is continually updated during 
the trimming procedure. It is stored under the array name RFE and is computed in 
subroutine BODYC. The sign convention for the inputs are given below noting that the 
in the CAMRAD.Modl inertial frame +X is in the forward direction while +z is down: 

AXTRIM = + forward = +a x 
AZTRIM = + up = -a ? 

Now, the accelerations in the body frame of reference may be written: 




AXTRIM 

"v/> 

► = [/?F£]< 

0 



AZTRIM 


or symbolically, 

= RFE{ 1,1) * AXTRIM - RFE( 1,3) * AZTRIM 
a yh = RFE{ 2,1) * AXTRIM - RFE{ 2,3) * AZTRIM 
a :h = RFE{ 3,1) * AXTRIM - RFE( 3,3) * AZTRIM 

The trim equations are written in terms of force rather than acceleration, therefore these 
terms must be converted to forces and then normalized. Converting to forces is 
accomplished by multiplying the accelerations by the vehicle mass. The conversion to 
forces and the normalization are both accomplished using the following conversion 
factor: 


F n = ma x = 
CF„ 


GW 


a. 


V g ) 

F.. 


r GW ^ 


o pnR 2 (FlR)~o 


\ g ) 


PkR 2 (CIR) 2 g 
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The normalization parameters are already available in CAMRAD.Modl under the 
symbolic names 


HMASS= 


( GWO 


V g ) 


DENSE= p 
FSCALE= Q 
RSCALE = R 
SSCALE= a 


The normalization factor in CAMRAD.Modl may now be written as: 

AFAC=HMASS/(DENSE*3.141593*RSCALE**2*(FSCALE*RSCALE)**2*SSCALE) 

The free flight trim equations in CAMRAD.Modl can now be modified using the force 
terms described above. The original force equations become (x-direction shown as an 
example): 



which is written in modified form as: 


( c \ 

^ FX 

(C FAX ) 

l O’ ) 

l"l C7 J 


or 


y [ Cpx ^ 


a J 


r 

I AX 


= 0 


where ^ contains the x-force contributions due to weight, rotor, fuselage, etc. 

C, 


<7 ) 


and 


a J 
accelerations. 


- contains the contributions to the x-force due to the specified inertial 


Appropriate printout has been added to the CAMRAD.Modl output to reflect the new 
capability. The modifications have been confined to the following subroutines: INITC, 
INPTN, PRNT, TRIMI, and TRIMP. The two common blocks TMDATA and TRIMCM 
have also been expanded to accommodate the additional data and these two appear in 
many subroutines. 

CAMRAD.Modl was modified to allow free flight trim option based on specified 
acceleration components. A preprocessor FORTRAN program determined the 
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acceleration components and these values were subsequently used in the 
CAMRAD.Modl input file. Pitch acceleration was ignored as described by the subtask. 
Three Crow’s Landing data sets were used for comparison with the numerical results. 
Two cases with 600 ft/min constant rate-of-descent and decelerations of 0.25 kt/sec and 
0.75 kt/sec were used. 

Case 1 Data Set 12563 600 FPM ROD W/0.25 KT/SEC DECEL FROM 90 TO 40 KIAS 

Characteristics of this data set obtained for a 600 ft/min rate-of-descent trajectory with a 
deceleration of 0.25 kt/sec are plotted in Figures 1 1 to 25. Figure 1 1 is a plot of the 
measured rate-of-descent. This data seems to indicate an average rate-of-descent of 800 
ft/min which is higher than the 600 intended. Figure 12 shows the actual trajectory of Z 
(ft) versus X (ft) and the chosen CAMRAD.Modl points that were selected for the 
comparisons are shown superimposed on this trajectory. This same information is shown 
as a (unction of time in Figure 13. The rate-of-descent can be calculated from this 
distance plot of the aircraft and it varies from 840 to 660 ft/min during the course of the 
trajectory. Additionally, as a function of airspeed, the descent from 64 kts (at which point 
the data is given) to 40 kts and lower is shown in Figure 14. The intended flight was 
supposed to start at 90 kts but this data was not in the data file for this data set. 
CAMRAD.Modl selected data points are superimposed on this graph at speeds 
distributed throughout the trajectory. Figure 15 shows the airspeed during the trajectory 
as a function ot time. Figure 16 shows the sideward flight of the aircraft during the 
trajectory. The magnitude is not that large over the course of 10000 ft indicating that the 
pilot was able to keep the aircraft on the course of the trajectory. The measured values of 
the angle-of-attack of the aircraft as measured from the aircraft boom is shown in Figure 
17 as a function of distance. The average value of 15° angle of attack seems to be too 
high for this flight. Hence, the angle of attack was not used in this comparison with 
CAMRAD.Modl values. This data was plotted with time in Figure 18 showing the same 
trend around 15 degrees. Roll was plotted as a function of distance in Figure 19. The 
value of roll was around 0° which is what the desired value should be as the aircraft trims 
to a zero roll value in actual flight. Figure 20 is the flight test data value of sideslip as a 
function of distance. It is very oscillatory in behavior as the pilot responds to the flight 
conditions encountered. Figure 21 shows the comparison of sideslip with time compared 
with results from CAMRAD.Modl. This comparison shows that the general trend of 
increasing sideslip with decrease in velocity with time is followed. At the larger velocity, 
the flight test required a higher sideslip value than that predicted. Sideslip is plotted 
with airspeed in Figure 22. As noted from Figure 21 in the plot with time, the 
comparison with sideslip at the higher speed was not matched for this flight test, probably 
owing to the local flight conditions at the time the flight test occurred which necessitated 
a higher sideslip. In general at the other speeds, the comparison with prediction is good. 
Figure 23 shows the aircraft angle of pitch with distance and the CAMRAD.Modl 
predictions. Although the aircraft pitch angle oscillates by about a degree in the flight 
test, the general trend of increasing pitch angle with lower speed as the pilot prepares to 
land is reflected by both sets of data. The comparison is generally good, considering that 
the flight test would have to compensate for any sudden gusts during the flight and these 
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are not simulated in the CAMRAD.Modl predictions. Figure 24 shows this same 
comparison with time and a similar trend is observed in this comparison as well. The 
oscillatory nature of the pitch amplitude is also present here as the pitch angle oscillates 
by one or two degrees during the flight trajectory. It is also interesting to look at the 
variation of pitch angle with airspeed as shown in Figure 25. The scatter in the pitch 
angle is apparent in this graph with a two degree variation at any speed indicating 
differing flight conditions during the trajectory as previously indicated. The 
CAMRAD.Modl predictions seem to fall in the general mean of the flight test data as 
would be expected by the ideal conditions used for the predictions. 

Case 2 Data Set 12565 600 FPM ROD W/0.75 KT/SEC DECEL FROM 90 TO 40 KIAS 

This case is the same constant rate-of-descent as the previous case 12563, with a larger 
deceleration rate of 0.75 kt/sec. Figure 26 shows the flight test data rate-of-descent 
measurement. Again it is around 800 ft/min on average throughout the trajectory with a 
large amplitude of oscillation of approximately 600 ft/min going to a maximum of 1400 
ft/min at some point. Figure 27 shows the trajectory in space with the CAMRAD.Modl 
points selected for this comparison shown on the graph. Figure 28 shows this same 
information as a function of time. These values indicate a rate-of-descent closer to 600 
ft/min than the previous case. The position with airspeed is shown in Figure 29. Actual 
airspeed during the trajectory is shown in Figure 30. In addition, it is interesting to note 
the variation in side motion given by the Y coordinate in Figure 31 shown with time. The 
actual flight trajectory confronts conditions which are not ideal, giving rise to some 
sideflight motion. Figure 32 shows the measured angle of attack during the trajectory. 

As in the previous case the average amplitude of the angle of attack of 15 degrees seems 
too large for this case and is therefore not being used during this comparison. The roll 
distribution with time is shown in Figure 33. The average value of roll here is a little 
over 1°, which is expected for a flight test trajectory when a 0 value of roll is difficult to 
maintain. This same distribution is shown with time in Figure 34. The distribution of 
roll with airspeed is shown in Figure 35. This distribution indicates a difficulty to 
maintain zero roll at the large starting velocities. It is possible that flight conditions at 
that time may have prevented trim to a value of roll close to zero. Sideslip for this 
trajectory is given in Figure 36 with CAMRAD.Modl predictions. For this case also, as 
in the previous case, there is a large oscillation in the sideslip value during the trajectory. 
Although the CAMRAD.Modl solution seems reasonable, gradually increasing in 
sideslip as the velocity decreases, the flight test data seems to indicate that the flight 
conditions were not as ideal as that assumed by CAMRAD.Modl. The same information 
is shown with time in Figure 37 and the same CAMRAD.Modl trend is seen here. 
Additionally, the sideslip variation with airspeed is shown in Figure 38 indicating that the 
agreement with prediction was good at the initial part of the trajectory and not so as the 
airspeed was decreasing. CAMRAD.Modl solutions and flight test comparison of 
aircraft pitch angle are shown in Figure 39. In this case, generally good agreement was 
obtained throughout the trajectory. A difference of approximately two degrees in pitch is 
acceptable for comparison with flight test under conditions which may not be ideal 
throughout the trajectory. Figure 40 shows this data as a function of time and similarly. 
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the agreement is fairly good. Pitch variation with airspeed is shown in Figure 41. In 
general, the comparison with CAMRAD.Modl results seems good. It appears that the 
comparison improves with lower velocity. 

Case 3 Data Set 12567 5.5 DEGREE APPROACH W/ 0.25 KT/SEC DECEL FROM 90 

TO 40 KIAS 

This case is a constant glideslope case of 5.5 degrees. Figure 42 shows the trajectory of 
height Z with distance X. It is approximately 5.5 degrees over the whole trajectory. 

Figure 43 shows the same information with time as the independent variable. Height as a 
function of airspeed is shown in Figure 44 showing the CAMRAD.Modl sample points 
for this trajectory. Figure 45 shows the airspeed variation of the trajectory with time. 

This trajectory shows some slight variations in the generally constant glideslope flight 
test. Figure 46 shows the side variation during the flight trajectory with time. During this 
flight there appears to be a large sideward flight component. This may have been the 
pilot compensating for flight conditions that are not ideal. Figure 47 shows the angle of 
attack variation with distance during this flight test and indicates an average of 15 degrees 
as in the previous two test cases chosen for this study. Again, this value seems too large 
for these flight conditions and therefore, angle of attack will not be used for this 
comparison. The roll angle with time is shown in Figure 48. For this case, the average 
roll is not 0 but approximately 2 degrees. In CAMRAD.Modl, all computations were 
carried out with the assumption of trim to 0 degrees of roll as in ideal flight conditions. 
This difference could contribute to any differences in flight test and computational 
results. In addition, the roll is seen to be quite oscillatory and attains a maximum value of 
approximately 6 degrees at an early point in trajectory. Figure 49 shows the variation of 
roll with airspeed. This graph more clearly shows the high value of roll at the initial large 
value of velocity in the trajectory. Figure 50 show the sideslip measured by flight test and 
that predicted by CAMRAD.Modl as a variation of distance. In general the trend of the 
CAMRAD.Modl results is correct in that the sideslip generally increases with decreasing 
velocity as the landing approaches. The flight test data is generally 5 degrees higher than 
predictions for this case. One can also observe the oscillations in the sideslip flight test 
data probably corresponding to varying flight conditions. Sideslip variation with time is 
shown in Figure 5 1 showing similar information to that just discussed. Figure 52 shows 
the variation of sideslip with airspeed during the trajectory. The general trend is indicated 
to be correct between the flight test data and the CAMRAD.Modl predictions. Figure 53 
shows the aircraft angle of pitch with distance from flight test data and CAMRAD.Modl 
solutions. The agreement here is good interrupted only by the large angle of pitch 
experienced by the aircraft during mid-trajectory. This is probably owing to flight 
conditions encountered during the constant glideslope trajectory. As the flight speed 
decreases toward the landing the pitch angle increases as expected for both sets of data 
presented. This information is presented as a function of time in Figure 54. The 
variation with airspeed can be seen in Figure 55. The agreement is good and shows the 
location of the CAMRAD.Modl points throughout the trajectory as a function of speed. 
Varying flight conditions are probably responsible for the scatter of data from the flight 
test. 
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Conclusion 


Analytical tools were developed to model noise abatement flight paths flown by 
Sikorsky’s S-76B helicopter. Sikorsky S-76B fuselage airloads data was included in 
CAMRAD.Modl. Rotor alone and full aircraft trim cases using CAMRAD.Modl were 
compared with wind tunnel and flight test data respectively and the comparison was 
satisfactory. A kinematic trajectory analysis tool was developed and interfaced with 
CAMRAD.Modl to be used in decelerating descents with constant glideslope and 
constant rate-of-descent trajectories. The comparison of the modified CAMRAD.Modl 
aircraft pitch and sideslip conditions generally agree with the Crow’s Landing flight test 
data. 
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FIGURE 6 MAIN ROTOR COLLECTIVE IN AIRCRAFT LEVEL FLIGHT 
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FIGURE 7 MAIN ROTOR LATERAL CYCLIC DURING LEVEL FLIGHT 
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FIGURE 8 MAIN ROTOR LONGITUDINAL CYCLIC DURING LEVEL FLIGHT 





FIGURE 9 FLOWCHART FOR APPROACH TRAJECTORY ACOUSTIC 
ANALYSIS 
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Equal Time Increments 
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FIGURE 10 CONSTANT GLIDESLOPE SPECIFIED BY EQUAL TIME 
INCREMENTS, CONSTANT RATE-OF-DESCENT SPECIFIED BY EQUAL 
RANGE INCREMENTS AND EQUAL TIME INCREMENTS 
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FIGURE 12 TRAJECTORY IN SPACE COORDINATES FOR TEST NUMBER 12563 (600 FT/MIN ROD) 
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FIGURE 13 TRAJECTORY IN TIME FOR TEST NUMBER 12563 (600 FT/MIN ROD) 
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FIGURE 14 TRAJECTORY WITH AIRSPEED FOR TEST NUMBER 12563 (600 FT/MIN ROD) 
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FIGURE 20 SIDESLIP DURING TRAJECTORY FOR TEST NUMBER 12563 (600 FT/MIN ROD) 
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FIGURE 23 AIRCRAFT PITCH ANGLE DURING TRAJECTORY FOR TEST NUMBER 12563 (600 FT/MIN ROD) 
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FIGURE 24 AIRCRAFT PITCH ANGLE (WITH TIME) DURING TRAJECTORY FOR TEST NUMBER 12S63 (600 FT/MIN ROD) 
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FIGURE 25 AIRCRAFT PITCH ANGLE (WITH AIRSPEED) DURING TRAJECTORY FOR TEST NUMBER 12563 (600 FT/MIN 

ROD) 
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FIGURE 28 TRAJECTORY IN TIME FOR TEST NUMBER 12565 (600 FT/MIN ROD) 
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FIGURE 30 AIRSPEED DURING TRAJECTORY FOR TEST NUMBER 12565 (600 FT/MIN ROD) 
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FIGURE 37 SIDESLIP TRAJECTORY FOR TEST NUMBER 12565 (600 FT/MIN ROD) 
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FIGURE 41 AIRCRAFT PITCH ANGLE (WITH AIRSPEED) DURING TRAJECTORY FOR TEST NUMBER 12565 (600 FT/MIN 

ROD) 
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FIGURE 43 TRAJECTORY IN TIME FOR TEST NUMBER 12567 (5.5 DEGREES APPROACH) 
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FIGURE 47 ANGLE OF ATTACK DURING TRAJECTORY FOR TEST NUMBER 12567 (5.5 DEGREES APPROACH) 
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FIGURE 48 ROLL DURING TRAJECTORY (FLIGHT TEST DATA) FOR TEST NUMBER 12567 (5.5 DEGREES APPROACH) 
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FIGURE 52 SIDESLIP DURING TRAJECTORY FOR TEST NUMBER 12567 (5.5 DEGREES APPROACH) 
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FIGURE 53 AIRCRAFT PITCH ANGLE DURING TRAJECTORY FOR TEST NUMBER 12567 (5.5 DEGREES APPROACH) 
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FIGURE 55 AIRCRAFT PITCH ANGLE (WITH AIRSPEED) DURING TRAJECTORY FOR TEST NUMBER 12567 (5 5 

DEGREES APPROACH) 





73 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data 
sources gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other 
aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and 
Reports. 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202*4302, and to the Office of Management and Budget, Paperwork Reduction Protect (0704-01 88) 
Washington, DC 20503 


1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 

February 2000 


4. TITLE AND SUBTITLE 

Rotoreraft Noise Abatement Flight Path Modeling 


6. AUTHOR(S) 

H. Murty and C. Berezin 


3. REPORT TYPE AND DATES COVERED 

Contractor Report 


5. FUNDING NUMBERS 

581-20-31-01 
NAS 1 -2(X)97 

Task 1 3 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Sikorsky Aircraft Corporation 
6900 Main Street 
Stratford. CT 06615-9129 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23681-2199 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA/CR-2000-209353 


11. SUPPLEMENTARY NOTES 

Technical Monitor: David A. Conner 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 

U nc lassi fled- U n 1 i mi ted 

Subject Category 71 Distribution: Nonstandard 

Availability: NASA CAS1 (301 ) 621-0390 


13. ABSTRACT (Maximum 200 words) 

This report addresses development of a rotor state/trim modeling capability for noise modeling of decelerating 
rotoreraft approaches. The resulting technique employs discretization of the descent trajectory as multiple 
steady state segments for input to CAMRAD.Mod 1 to predict rotor states for acoustic analysis. Deceleration is 
included by modifying the CAMRAD.Mod I tree flight trim options to allow trim to the specified acceleration/ 
deceleration components. 


14. SUBJECT TERMS 

Rotoreraft, Helicopter, Noise, Abatement, Acoustics, Modeling, Flight Path 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


15. NUMBER OF PAGES 

77 


16. PRICE CODE 

A05 


20. LIMITATION 
OF ABSTRACT 

UL 


NSN 7540-01-280-5500 


74 


Standard Form 298 (Rev. 2-69) 

Prescribed by ANSI Std. Z-39-18 
298-102 















